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Using sum rules, the dipolar terms can be eliminated from the commonly-used sum-over-states 
(SOS) expression for nonlinear susceptibilities. This new dipole-free expression is more compact, 
converges to the same results as the common SOS equation, and is more appropriate for analyzing 
certain systems such as octupolar molecules. The dipole-free theory can be used as a tool for 
analyzing the uncertainties in quantum calculations of susceptibilities, can be applied to a broader 
set of quantum systems in the three-level model where the standard SOS expression fails, and more 
naturally leads to fundamental limits of the nonlinear susceptibilities. 
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I. INTRODUCTION 

The sum-over-states (SOS) expression for the 
nonlinear-optical susceptibilities, which are expressed 
in terms of the matrix elements of the dipole operator, 
—ex nm , and the energy eigenvalues, E n , have been 
used extensively over the last 4 decades as a theoretical 
model of the nonlinear response as well as a tool for 
analyzing experimental dispersion data. Indeed, the 
two-level model has guided the development of molecules 
with large second-order nonlinear-optical susceptibilities 
(also called hyperpolarizabilities) . Similarly, the SOS 
three-level model for the third-order susceptibility (also 
called the second hyperpolarizability) has led to an un- 
derstanding of the nature of the states, and interactions 
between them, that yield the largest response. 

While the SOS expression has reigned supreme for 4 
decades, there are several critical issues that have never 
been addressed. Because the expression is over-specified, 
redundant information is required to calculate nonlinear 
susceptibilities. This redundancy not only leads to inef- 
ficiencies in the computational process; but, the SOS ex- 
pression can lead to reasonable-looking results even when 
unphysical parameters are used as the input. As such, 
the underlying physics of the nonlinear-optical response 
may be misinterpreted, leading to erroneous conclusions. 

The sum rules are quantum mechanical identities that 
relate the dipole matrix elements and energies to each 
other; so, the SOS hyperpolarizability can be expressed 
in terms of a subset of the dipole matrix. Q In this work 
we show that all the dipolar terms can be eliminated to 
yield a simplified expression that is equivalent to the full 
SOS expression. This theoretical result can be used as a 
tool for studying the nonlinear-optical response. For ex- 
ample, differences between the full SOS expression and 
the dipole-free expression can be used to estimate the 
uncertainties in quantum calculations since such discrep- 
ancies are a sign that the sum rules have been violated 
(due to incorrect dipole matrix elements and energies, or 
truncation errors). The approximate solution of a par- 



ticle in a tilted box is analyzed in this way as an illus- 
tration. More importantly, the three-level model of the 
dipole-free theory may be applicable when the standard 
three-level SOS expression fails, thus providing a theoret- 
ical tool that covers a broader base of quantum systems. 
To test our new theoretical expression we show that for 
a clipped harmonic oscillator, for which analytical solu- 
tions to the Schrodinger Equation can be computed, the 
two theories converge to the same results in the limit of 
an infinite-level model. Finally, we shall show that the 
dipole-free expression more naturally leads to fundamen- 
tal limits of the nonlinear susceptibilities. 



II. THEORY 

In this section, we show how the SOS expression of the 
first and second hyperpolarizabilities are simplified using 
the sum rules. The commonly-used SOS expression for 
any diagonal component of (3 is given 
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where e is the magnitude of the electron charge, x nm the 
n, m matrix element of the position operator, Ax 71 q = 
Xnn — %oo is the difference in the expectation value of 
the electron position between state n and the ground 
state 0, E n o = E n — Eq is the energy difference between 
the excited state n and the ground state, E m a(tj}p) 
/'.*,,; i — tkop, and Tiujp is the energy of one of the photons. 
The primes indicate that the ground state is excluded 
from the sum and the permutation operator P(ui a ,uip) di- 
rects us to sum over all six frequency permutations given 
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by the Feynman Diagrams. Since the dipole moment of 
the molecule is proportional to the position (p x = —ex), 
we loosely call x nm the transition moment and x nn the 
excited state dipole moment. The first term in Equation 
13 is called the dipole term and the second term the oc- 
tupolar term; and, as we shall see below, the dipole term 
can be expressed in terms of the octupolar one using the 
sum rules. 

The generalized Thomas-Kuhn sum rules are a direct 
consequence of the Schrodinger Equation (without any 
approximations) and relate the matrix elements and en- 
ergies to each other according to: [j| 
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where m is the mass of the electron, and N the number 
of electrons in the molecule. The sum, indexed by n, is 
over all states of the system. Equation [3 represents an 
infinite number of equations, one for each value of m and 
p. As such, we refer to a particular equation using the 
notation (m,p). 

To eliminate the dipole term, we consider the Equation 
(jn,p) with m 7^ p: 
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Equation 13 can be rewritten by explicitly expressing the 
n = m and n = p terms: 
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Setting p = in Equation[5]and solving for Ax n o \xon \ 
after multiplying through by XQ m , we get 
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Substituting Equation with m <-> ro into Equation [21 
we get the final result, 
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We call this form of /3 the dipole-free expression or the 
reduced hyperpolarizability. The second term in brack- 
ets is the dispersion term that results when the dipolar 
terms are eliminated. In the standard SOS expression, 
the simplest approximation is the two-level model, with 
parameters xio, Axio, and E\$. The simplest approxi- 
mation to Equation [5] is the three-level model with pa- 
rameters xio, X20, x\2, E10, and i?2o- This is in contrast 
to the standard SOS expression, where the three-level 
model has two additional dipole terms. 

It is important to note that while the dipole-free ex- 
pression may seem to be less general than the common 
SOS one, when all states arc included, it is fully equiv- 
alent. Because the sum rules are a direct consequence 
of the Schrodinger Equation, they can not be violated in 
any system, be it an atom, molecule, or crystal. The SOS 
expression, in both forms, can be evaluated for unphys- 
ical values of the matrix elements - yielding nonsensical 
values of the hyperpolarizability. However, in its nonre- 
stricted form, there is more room for introducing errors. 
The restriction imposed on the SOS expression used to 
get the dipole-free equation prevents certain unphysical 
combinations of dipole and octupolar terms, so in non- 
truncated form, is more robust. For example, the lowest 
truncated-state model in the standard SOS expression 
is the two-level model, which only describes transitions 
in which the dipole moment changes between these two 
states. However, it ignores all octupolar terms. On the 
other hand, the dipole-free expression - when truncated 
- can approximate molecules with octupolar character 
[IS El as measured with hyper- Rayleigh scattering; 0,0 
as well as dipolar terms, which are implicitly taken into 
account by the extra dispersion term in the reduced hy- 
perpolarizability, as given by Equation [S] 

All higher-order nonlinear susceptibilities can be 
treated in the same way. As an illustration, we briefly 
consider the third-order susceptibility. For any diagonal 
component of 7, the second hyperpolarizability, along the 
x-direction is given by: 
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where Di(n, m, ...) are energy denominators and n, m, ... 
are arguments that show which energies are represented 
(i.e. E n o, E m Q, ...). In analogy to (3 as given by Equa- 
tion [21 the denominators are of the form Di(n,l,m) — 
fi 3 'E n o(u} a )E m o(ui/3)Eio(u)s) /4e 4 and are most easily de- 
termined using Feynman Diagrams for the particular 
phenomena of interest. There are two terms in Equa- 
tion IS] that depend on the dipole moment, which can be 
expressed as, 
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Using Equation with p = and n = I, Equation El 
becomes 
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Using Eg uat ions 1111 and 1131 Equation[S]for 7 can be writ- 
ten in dipole-free form in analogy to Equation [SJ Given 
the algebraic messiness of the result, it is not presented 
here. 



III. APPLICATIONS 

In this section, we apply the dipole-free theory to sev- 
eral problems to show both its usefulness; and, to confirm 
that in the limit of including all states in the SOS expres- 
sion, our new theory and the standard SOS expression 
converge. Furthermore, we show that when approximate 
wavefunctions are used, such as the particle in a tilted 
box, the two results do not converge, showing how such 
a comparison can be used to assess the accuracy of the 
calculated nonlinearities. In addition, we show how the 
reduced hyperpolarizability leads to a more elegant cal- 
culation of the fundamental limits of nonlinear suscepti- 
bilities. 

First, we apply the theory to the calculation of the fun- 
damental limits of the hyperpolarizability. 0- El fill 
12] The issue of fundamental limits has been an impor- 
tant one since it guides the applied researcher in making 
better materials and devices while giving the theorist a 
method of understanding the nuances of what makes a 
large nonlinear-optical response. El Using the new 
theory, we show that the results are the same; but, leads 
to a more elegant approach that illustrates the equiva- 
lence of viewing a molecule in terms of the standard ex- 
pression that includes dipole terms or the reduced form 
in terms of octupolar terms. This provides an important 
paradigm shift in the sense that the sum rules show that 
the two limiting cases are closely related while the general 
nonlinear-optics community operates on the assumption 
that the two are independent. 

Second, the clipped harmonic SOS expression - for 
which exact analytical wavefunctions can be calculated - 
the SOS expression and the dipole-free theory presented 



here are compared. The fact that the two converge shows 
that the two expressions are identical in the infinite-level 
model. It is common for theorists and experimentalists 
to use limited-state models. The two- and three-level 
models have been successfully applied to understanding 
the dispersion and magnitude of the second and third- 
order susceptibilities, but clearly, such a simplified view 
can not be universally correct. Indeed, the three level 
model for the SOS expression and the dipole free one are 
totally different functions with different dispersion. 

There are two important ramifications of this observa- 
tion. First, in cases where the standard truncated SOS 
expression is inconstant with observation, the dipole-free 
expression may be more appropriate. As such, the dipole- 
free expression provides researchers with a tool to study 
a class of systems that was previously unaccessible. Oc- 
tupolar molecules may be one such class. Secondly, a 
comparison of the dipole-free expression with the stan- 
dard SOS theory can be used to estimate whether or not 
a theoretical calculation has converged without the need 
for including more and more numbers of states to test 
for convergence. So, the theory presented here can save 
on computational time while providing the theorist with 
another tool. 

Third, /3 of a particle in a tilted box is analyzed. Since 
this system is solved with perturbation theory, it is pos- 
sible to study how small errors in the matrix elements 
affect the nonlinear response predicted by the two the- 
ories. Most quantum chemical calculations yield only 
approximate wavefunctions, so this example shows how 
the two expressions can yield different results even in 
an infinite state model. It must be stressed that there 
are no means for determining which calculation yields 
the "correct" (3 values that would be found experimen- 
tally. Perhaps more importantly, a comparison between 
the two theories provides an estimate of the uncertainty 
of the calculation: The calculated values of the nonlin- 
earity can only be trusted to within the range between 
the two values. If both converge to greatly differing val- 
ues in the infinite-state limit, then this suggests that the 
wavefunctions may be unphysical. 



A. Fundamental Limits 

To illustrate the usefulness of the dipole-free SOS ex- 
pression, we apply it to calculate the fundamental upper 

limit of /3HHlll!3,[ll[I2- We start with the sum rulcs 
(0, 0) and (1, 1) truncated to three levels, which yield 
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respectively, where 



„MAX 1 2 



2mE%o 



-N. 



(15) 



Substituting Eg uat ions 1131 and 1 1 41 into Equation [8] in the 
off- resonance limit (uj a = u)p = Q), we get 
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where X = x w /x^ AX and E = E w /E w . 

G and X are maximum at G( ~\/3) = 1 and /(0) = 1, 
yielding, 
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This is identical to the results from the usual sum-over- 
states expression; however, the calculation is much more 
concise and elegant because the dipolar term does not 
need to be considered. 



B. The Clipped Harmonic Oscillator 

In this section, we test the dipole-free SOS expression 
by comparing the results it gives with the standard SOS 
expression for a potential in which the Schrodinger Equa- 
tion can be solved analytically. This approach ensures 
that the energies and dipole matrix elements are physi- 
cally sound and that pathologies or inaccuracies normally 
inherent in approximation techniques are avoided. We 
use the exact solution to the clipped harmonic oscilla- 
tor (CHO) (where the potential is harmonic for x > 
and infinite for x < 0) since it is the simplest case of an 
asymmetric potential that yields a large hyperpolariz- 
ability that is in fact near the fundamental limit. |15j| The 
matrix elements of the position operator of the clipped 
harmonic oscillator (CHO) are given by, 
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where the dimensionless matrix g mn is defined by 
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FIG. 1: /3(0)//3q /ax , the zero- frequency (off-resonance) limit 
of ft - normalized to the off-resonant fundamental limit - as a 
function of the number of excited states included in a clipped 
harmonic oscillator and ID tilted Box for the standard SOS 
model and the dipole-free SOS expression. 



where to and n are restricted to the odd integers. The 
energy of state n is given by 
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Figure ^ shows the calculated off-resonant hyperpo- 
larizability normalized to the maximum off-resonant hy- 
perpolarizability as a function of the number of states 
included in the calculation. Both theories converge to 
the same result as the number of states included in the 
sums is large, showing that the two models are identi- 
cal. Note that the standard SOS expression converges 
more quickly than the dipole-free expression, which sug- 
gests that the clipped harmonic oscillator is more dipo- 
lar in nature. Presumedly, an octupolar molecule would 
be better modelled with the dipole-free term, resulting 
in faster convergence; though, there are no simple ex- 
actly soluble octupolar potentials. The average of the 
two models is also shown, suggesting that a variational 
principle applied to a weighted average (with the weights 
as parameters) may yield the exact result with only a few 
terms. 

Figure |21 shows the dispersion predicted by both mod- 
els for a CHO in the 3- and 6-level models for the second 
harmonic generation hyperpolarizability as a function of 
the energy of the fundamental photon. The two the- 
ories agree well in the 6-level model except near reso- 
nance. The insets show an expanded view of the regions 
in which the two theories disagree the most. In the 25- 
level model(Figure I3J , the agreement is excellent at all 
wavelengths, as expected since the CHO is an exact so- 
lution to the Schrodinger equation. 
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FIG. 2: The normalized second harmonic hyperpolarizability 
((3(E)/(3q IAX ) as a function of the incident photon energy for 
a 3- and 6-level model of a clipped harmonic oscillator for 
standard and dipole-free SOS expressions. Insets show mag- 
nified view of key regions as indicated by the dashed boxes. 
The first excited state energy is arbitrarily set to leV. 
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FIG. 4: The normalized second harmonic hyperpolarizability 
(/3(E)//3q IAX ) as a function of the incident photon energy 
for a 25-level model of a particle in an asymmetric box for 
standard and dipole-free SOS expressions. The inset shows a 
magnified view. 
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FIG. 3: The normalized second harmonic hyperpolarizability 
(/3(E)/(5q IAX ) as a function of the incident photon energy for 
a 25-level model of a clipped harmonic oscillator for standard 
and dipole-free SOS expressions. The inset shows a magnified 
view. 



C. Particle in a Tilted Box 

Next we consider a particle in a 1-dimensional box that 
is perturbed by the potential V = ex to make the system 
asymmetric. First-order perturbation theory is used to 
get the wavefunction to first-order in e, from which the 
matrix elements of x are calculated. (3 is calculated from 
these matrix elements also to first-order in e. This is 
an interesting example because the wavefunctions, while 
reasonably accurate, are nevertheless only approximate. 
Figure^shows /?/ (3q IAX for the two models as a function 
of the number of states and Figure 01 shows the 25-level 
model. Note that the matrix elements are accurate to 



better than 5%, yielding convergence of the off-resonance 
limit of the two 25-level models of j3 to better than 7% 
of each other. However, near resonance, the two models 
do not agree as well quantitatively - though the quali- 
tative features are similar. These variations are due the 
inaccuracies introduced by the approximations used in 
calculating the wavefunctions, so it is not possible to de- 
termine which model is more accurate. However, based 
on the two dispersion graphs, it is clear that the dipole 
free-expression and standard SOS expressions are equiv- 
alent to within the levels of uncertainty one expects from 
the level of approximation used. 



IV. CONCLUSION 

In conclusion, we have derived an expression that is 
equivalent to the standard SOS equation for [3 and 7, 
but does not include dipole terms. The fact that they are 
identical is illustrated with the exact wavefunctions of a 
clipped harmonic oscillator; when the number of terms 
included in the sums is large, the two results converge. 
However, when the approximate wavefunctions of a par- 
ticle in a tilted box are used, the two expressions do not 
converge, illustrating how the difference can be used to 
estimate the uncertainty in the result. Furthermore, such 
a variance may also be a sign that the wavefunctions used 
violate the sum rules. 

The dipole-free expression is more compact; and, when 
truncated to a finite number of states is easier to apply 
to certain classes of problems, such as calculating the 
fundamental limits of the nonlinear susceptibility. The 
dipole-free expression is a new tool for studying classes 
of molecules that are not well described by the truncated 
SOS expression. As such, in may, for example, provide a 
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more accurate means for analyzing the dispersion of [3 for 
octupolar molecules. The standard approach is to trun- 
cate the sums in EquationElto the first two excited states 
(yielding the term with numerator x$\Xi2X2o)- Clearly, 
the dipole-free 3-level model includes more information 
so may be a more accurate expression for the dispersion 
than simply setting the dipole term in Equation|21to zero 
in the standard 3-level model. 

The new theory for /3 and 7 presented here is there- 
fore an additional avenue for analyzing molecules that go 
beyond the common dipolar push-pull paradigm, can be 



used to assess the accuracy of molecular orbital calcula- 
tions, and sheds new light on the fundamental limits of 
nonlinear susceptibilities. 
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